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THE CONNECTION BETWEEN JETS, ACCRETION DISKS, AND BLACK HOLE MASS IN BLAZARS 



^1" 

o 
o 

(N 
O 

Q 



> 
o 

o 



Of 
i 1 

o 



X 



Jian-Min Wang , Bin Luo and Luis C. Ho 

Received 2004 August 1 7; accepted 2004 September 22; published October 4 

ABSTRACT 

The properties of relativistic radio jets are thought to be closely connected with the properties of accretion disks 
in active galactic nuclei. We explore this issue using a sample of 35 blazars with very-long baseline observations, 
for which we can estimate the kinetic powers of their relativistic jets, the bolometric luminosity of their accretion 
disks, and the masses of their black holes. Contrary to previous claims, we find that the jet kinetic power is 
significantly correlated with the disk luminosity. This supports the notion that the disk is somehow coupled to 
the jet, even on parsec scales. Moreover, we show that the correlation improves by including the black hole mass 
as a second parameter. The dominance of the jet, parameterized as the ratio of the jet kinetic power to the disk 
luminosity, is largely controlled by and is inversely correlated with the Eddington ratio of the accretion disk. This 
empirical relation should serve as a useful guide for theoretical models for jet formation. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — (galaxies:) BL Lacertae 
objects: general — galaxies: nuclei — (galaxies:) quasars: general 



estimated from modeling of their spectral energy distribution. 
They found that the jet power is linearly proportional to the disk 
power, with flat-spectrum radio quasars radiating at a greater ef- 
ficiency than BL Lac objects. In view of the importance of this 
issue, it certainly merits reexamination. 

Jet formation may depend on the state of the accretion disk. 
Marscher et al. (2002) found that the radio galaxy 3C 120, as 
in microquasars, ejected superluminous radio blobs following 
a dip in the X-ray emission. It is essential to probe the accre- 
tion rates of the disks to test the jet-disk connection. Wang et 
al. (2002, 2003, 2004) showed that most BL Lac objects have 
optically thin ADAFs, while most radio-loud quasars have stan- 
dard accretion disk. It is interesting to note that the degree of 
radio-loudness in low-luminosity active galactic nuclei (AGNs) 
is determined by the Eddington ratio (Ho 2002), and BH accret- 
ing systems obey a "fundamental plane" defined by their radio 
emission, hard X-ray emission, and BH mass (Merloni, Heinz, 
& Di Matteo 2003). The highly anisotropic emission from the 
jet depends on the dissipation of the jet kinetic energy, and the 
radiation efficiency may be controlled by both the state of the 
central engine and the environment of the jet (Celotti & Fabian 
1993; Wang, Staubert, & Ho 2002). However, how the kinetic 
luminosity, a more fundamental parameter of the relativistic jet, 
depends on the properties of the accretion disk is not yet known. 

In this paper, we tackle this issue with the sample of blazars 
from CPG. We show that in blazars there is a well-defined cor- 
relation between kinetic luminosity, BLR luminosity, which we 
use as a surrogate to represent the bolometric luminosity of the 
accretion disk, and the BH mass. Moreover, the degree of jet- 
dominance diminishes with increasing Eddington ratio. 

2. THE SAMPLE 

For the purposes of the present paper, we need the kinetic 
luminosity, bolometric luminosity, and the BH mass. Ac- 
cording to CPG, the kinetic luminosity is given by L^ m = 
fTrRy Lm n e m e c 2 T 2 j3c, where T^vlbi is the size of the radio jet, 
p is the jet velocity in units of the light speed (c), n e is the num- 
ber density of relativistic electrons, T is the Lorentz factor of 
the jet, and m e is the electron mass. 7?vlbi is an observable pa- 
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1. INTRODUCTION 

Many of the main properties of blazars can be explained by 
relativistic jets (see review by Urry & Padovani 1995); how- 
ever, jet formation remains one of the unsolved fundamen- 
tal problems in astrophysics (e.g., Meier, Koide, & Uchida 
2001). Three scenarios have been suggested for the origin of 
jets: (1) from the vicinity of the black hole (BH) event hori- 
zon via extraction of the spin energy of the BH (Blandford & 
Znajek 1977); (2) from a wind from a magnetized accretion 
disk (Blandford & Payne 1982); and (3) from an advection- 
dominated accretion flow (ADAF; Narayan & Yi 1994), either 
dominated by outflows (Blandford & Begelman 1999) or a jet 
(Falcke, Kording, & Markoff 2003). Though there are still 
some unresolved issues with each of these models, it is clear 
that an accretion disk is involved in all of them. It is expected 
that there is a connection between the relativistic jet and the 
accretion disk. 

The jet-disk connection has been extensively explored by 
many authors and in different ways. Since emission from the 
relativistic jets is highly beamed, one needs to use isotropic ra- 
diation to investigate the connection. On large scales (^kpc- 
Mpc), many authors have tested the correlation between differ- 
ent emission-line luminosity and radio luminosity (Stockton & 
MacKenty 1987; Baum & Heckman 1989; Rawlings & Saun- 
ders 1991; Falcke, Malkan, & Biermann 1995; Xu, Livio, & 
Baum 1999; Cao & Jiang 2001; Wang, Ho, & Staubert 2003). 
Indeed, a strong link between the jet and the disk has been 
found, although we should note that it is indirect because the ex- 
tended radio emission and the line emission have quite different 
lifetimes. On ~pc scales, radio data from very long-baseline 
interferometry (VLBI) have been used to calculate the kinetic 
luminosity to test the jet-disk connection (Celotti, Padovani, & 
Ghisellini 1997; hereafter CPG). CPG found no correlation be- 
tween the broad-line region (BLR) luminosity (Lblr) and the 
kinetic luminosity (LkinX although there are some suggestive 
hints of a jet-disk connection. Maraschi & Tavecchio (2003; 
hereafter MT03; see also Maraschi 2001) investigated the prob- 
lem using a sample of blazars for which the jet power could be 
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REFERENCES/NOTES:- C92: Corbin 1992; CPG: Celotti et al. 1997; F03: Falomo, Carangelo, & Treves 2003; L96: Lawrence et al. 1996; M96: Marziani 
et al. 1996; MC04: Marchesini & Celotti 2004; 084: Oke, Shields, & Korycansky 1984; SF97: Scarpa & Falomo 1997; SFK89: Stickel, Fried, & Kiihr 
1989; SKF93: Stickel, Kiihr, & Fired 1993; SFK93: Stickel, Fried, & Kiihr 1993; SK93: Stickel & Kiihr 1993; SM:87 Stockton & MacKenty 1987; SS91: 
Steidel & Sargent 1991; UOO: Urry et al. 2000; V95: Vermeulen et al. 1995; W86: Wilkes 1986; W95: Wills et al. 1995; WU02: Woo & Urry 2002. 



HPQ: core-dominated high-polarization quasars; LPQ: core-dominated low-polarization quasars; LDQ: lobe-dominated quasars; BL: BL Lacs. 
°A11 the data are scaled to luminosity distances calculated assuming Hq = 75 km s~' Mpc~', Q m = 0.3, and Q\ = 0.7. *Doppler factor. c Mrk 421 has line 
flux for Ha instead of H/3 or Mg II . H/3 or Mg II lines. The upper limit on its core size given in Ghisellini et al. (1993) gives only an upper limit on logLyn. 
The value of logLki n listed is based on the new VLBI observations of Piner et al. (1999). rf Mrk 501 has line flux for Ha instead of H/3 or Mg II . 



rameter, while n e and T can be obtained by fitting the observed 
spectrum. However, the largest source of uncertainty for the ki- 
netic luminosity lies in the factor /, which describes the particle 
composition of the jet. For a jet composed of a proton-electron 
plasma, / has a value of ~ 1840, while for a pair plasma / = 1. 
As in CPG, we take / = 1 . In the present paper, we use a sub- 
set of the sample of CPG, which provides Lblr and Ly n . Our 
subset is limited by the availability of the following data: (1) 
£kin, (2) BH mass estimate, and (3) flux for the H/3 A4861 or 
Mg II A2800 line, which, as described below, we use to esti- 
mate Lblr- We omit those sources with an upper limit on Ly n 
because they only have an upper limit on /?vlbi- Ghisellini et al. 
(1993) give the Doppler factors of the sample, estimated from 
the classical limits of synchrotron self-Compton emission. In 
this study, we estimate BH masses using one of three methods, 
depending on the availability of the necessary data: (1) the rela- 
tion between stellar velocity dispersion and BH mass, as fitted 
by Tremaine et al. (2002); (2) the linewidth-luminosity-mass 
scaling relation, which is calibrated against local AGNs stud- 
ied using reverberation mapping (Kaspi et al. 2000); and (3) 
the correlation between host galaxy luminosity and BH mass, 



as given by McLure & Dunlop (2001). When the linewidth 
(FWHM) of H/3 is not available, we use the FWHM of C IV or 
Mg II as suggested by Vestergaard (2002) and McLure & Jarvis 
(2002), respectively, to estimate the BH mass. Our paper ana- 
lyzes 35 blazars (Table 1) from the sample of CPG. When the 
BH mass is estimated from the linewidth-luminosity-mass rela- 
tion, we use an extrapolation of the continuum from the magni- 
tudes given in Veron-Cetty & Veron (2003), assuming a spec- 
trum of the form F v oc v" 5 . 

The uncertainty in the M B h - c relation is small, < 0.21 dex 
(Tremaine et al. 2002). On the other hand, the uncertainty 
on the zeropoint of the linewidth-luminosity-mass relation is 
approximately 0.5 dex (Gebhardt et al. 2000; Ferrarese et al. 
2001). We note, however, that this method may systematically 
overestimate the BH masses in radio-loud quasars if the optical- 
ultraviolet continuum in these objects is significantly beamed. 
Unfortunately, it is difficult to estimate the magnitude of this 
effect, since the amount of beaming is wavelength dependent 
(e.g., Bregman 1990), and we have no good constraints for 
the optical-ultraviolet region. According to McLure & Dun- 
lop (2001), the M B h -M R correlation for quasar host galaxies 
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FIG. 1. — Correlation between jet kinetic luminosity and (a) Lblr and (b) Lblr and Lgdd for BL Lac objects (stars), lobe-dominated quasars (open circles), 
core-dominated, low-polarization quasars (squares), and core-dominated, high-polarization quasars (filled circles). The uncertainties on Lblr are based on the 
standard deviation of Mg II /La and H/3/La in the composite quasar spectrum of Francis et al. (1991), while the uncertainties on L^in are assumed to be 0.5 dex for 
all sources. The error bars for the abscissa in (b) take into account the uncertainties in Lblr and Mbh (see §2). 



has an uncertainty of 0.6 dex. Given the intrinsic uncertainty of 
the different BH mass estimators and the heterogeneity of our 
sample, we estimate that the individual BH masses may have 
an uncertainty as large as ^1 dex. 

3. EMPIRICAL CORRELATIONS 

The BLR luminosities given in CPG were derived by scaling 
several strong emission lines to the quasar template spectrum 
of Francis et al. (1991), using Lya as a reference. With this ap- 
proach, CPG find no significant correlation between Lblr and 
Lkin- We note, however, that some of the key emission lines 
used in their method have a very large scatter in their ratio with 
respect to Lya. For example, C IV A1550/Lya = 63 ±41 and 
Mg II A2800/Lya = 34 ± 20. The variation in the strength of 
C IV and Mg II relative to Lya stems from genuine object-to- 
object variations in the detailed physical conditions of the BLR 
in quasars. Since this scatter propagates directly into an error 
in Lblr, it is possible that this is responsible for the apparent 
lack of a correlation between Lblr and Ldn- On the other hand, 
the ratio H/3/Lya is largely governed by recombination physics 
(Netzer 1990) and hence is much better defined. Francis et al. 
(1991) find H/3/Lya = 22 ±4.1. Here, we will adopt H/3 alone 
as the fiducial line to estimate Lblr; when H/3 is not available, 
we use Mg II , which has the second smallest scatter relative to 
Lya. Following the notation of CPG (their Eq. 1), 

(1) 



Lblr : 



16.35LMg ii, 



Les,(MgII) 

where (L* LR ) = 555.77, L est (H/3) = 22, and L est (MgII) = 34. Ta- 
ble 1 gives the values of Lblr estimated from Equation 1 . 

Figure la shows a strong correlation between L B lr and L kin . 
We have tested the partial correlation among Ld n , Lblr, and 
c/ L , the luminosity distances of the objects, so as to avoid the 
possible spurious effects caused by their common dependence 
on the luminosity distance. Using the least square method of 
multivariate regression, we find that logLdn = 0.431ogLBLR + 
O.lllogc/L, with partial coefficient r^ L = 0.4 and r£ BLR = 0.62. 
This test shows that the correlation between Ld n and L B lr is 
intrinsic rather than introduced by the dependence of <f L . 



Since we are exploring the possible dependence of jets on 
the accretion disks (using Lblr as a substitute for the disk bolo- 
metric luminosity), L B lr should be treated as the independent 
variable. The least square method gives the following regres- 
sion fit in Figure la (excluding the prominent outlier 3C 345): 

logL km = (0.37 ±0.08)logL BLR + (29.78 ±3.7), (2) 

with Pearson's coefficient r = 0.62 and a probability of p = 
1.0 x 10~ 4 for rejecting the null hypothesis of no correlation. 
Assuming that the BLR luminosity is due to the reprocessing of 
the primary continuum radiation from the accretion disk, Equa- 
tion 2 can be interpreted as evidence for a strong connection 
between the relativistic jet and the accretion disk. The relation 
L^ oc L B lr is much flatter than the linear relation between jet 
power and disk power found by MT03. It is unclear what the 
source of discrepancy is. We simply note that our procedure for 
estimating disk luminosities differs from that of MT03. MT03 
obtain their disk luminosities either directly from the optical- 
UV luminosity of the big blue bump or from the original pre- 
scription of CPG. 

Despite the strong correlation exhibited in Figure la, clearly 
there is still significant scatter. The scatter around the best-fit 
line is a = 0.37 dex. Although a significant fraction of this scat- 
ter may be due to measurement and/or systematic errors in Ld n 
and Lblr, it is possible that it indicates that a second parame- 
ter is important. A nature choice to investigate is the BH mass, 
here cast in terms of its Eddington luminosity, L E( jd. Using the 
least square method of multivariate regression, we indeed find 
that Lyn correlates with both L B lr and L E dd- Excluding 3C 345 
and 0235+164 (its BH mass is an upper limit), we obtain 



logLkin = (0.25 ± 0.09) logL BL R+ 

(0.65 ± 0.25) logL Edd + 5 .07 ± 10.05 , 



(3) 



with Pearson's coefficient r = 0.70. The new correlation (Fig. 
lb) shows a smaller scatter, which has reduced to a = 0.31 dex. 
The signifcance of adding the Eddington luminosity as a new 
parameter in the correlation (3) can be tested by the L-test de- 
fined by Bevington (1969). We find F x = 7.36 and probability 
p = 1 .09 x 10~ 2 . The improvement in the correlation by includ- 
ing L Ea d as an additional variable suggests that the jet power 
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depends on both the disk luminosity (or accretion rate) and BH 
mass. Whether additional variables are required, however, can- 
not be decided on the basis of the data presented here. The 
scatter in Figure lb, though reduced, remains substantial, but 
the significant uncertainties currently associated with the ob- 
servables prevent us from assessing how much of the scatter is 
intrinsic. 

Equation 3 can be casted in a different form. Defining the 
"jet-dominance" factor as J 7 ] = Ljdn/^boi, the Eddington ratio as 
£ = Lboi/^Edd, and assuming Lboi ~ IOLblr (Netzer 1990), it 
follows that 

logJj = -0.101ogL bol -0.651og£ + 4.82. (4) 
This implies that the relative importance of the jet power com- 
pared to the disk luminosity is mainly controlled by, and is 
inversely dependent on, the Eddington ratio. MT03 reached 
a qualitatively similar conclusion in their analysis. This rela- 
tion provides a physical basis for understanding the spectral se- 
quence of blazars (Fossati et al. 1998; Wang et al. 2002, 2004; 
MT03). This empirical relation may also help guide theoretical 
models on jet formation. 

4. CONCLUSIONS 

Using rough estimates of the BH masses and bolometric lu- 
minosities for a sample of 35 blazars with VLBI observations, 



we find that the kinetic power of the relativistic jet in these ob- 
jects is strongly coupled both to the accretion disk luminosity 
and the BH mass. Equivalently, the degree of jet dominance, 
defined as the ratio of the jet kinetic power to the bolometric 
luminosity of the disk, is largely controlled by and is inversely 
correlated with the Eddington ratio of the accretion disk. This 
empirical relation should provide strong constraints on theoret- 
ical models for jet formation and jet-disk coupling. 

There is growing evidence in recent years that BH accre- 
tion in AGNs provides an important source of energy feedback 
to their host galaxies and their large-scale environments (see 
review by Begelman 2004). Relativistic jets from radio-loud 
sources, in particular, may be particularly effective in heating 
the intracluster medium in galaxy clusters. The empirical cor- 
relations described in this study, if they can be properly cali- 
brated, opens up the possibility of using them as tools to inves- 
tigate quantitatively accretion-powered feedback in the context 
of galaxy formation. 
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